Rossby wave equilibria and zonal jets 
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The problem of coherent vortex and zonal jet formation in a system of nonlinear Rossby waves 
is considered from the point of view of the late time steady state achieved by free decay of a given 
initial state. Statistical equilibrium equations respecting all conservation laws are constructed, 
generalizing those derived previously for 2D inviscid Euler flow. Jet-like solutions are ubiquitous, 
with large coherent vortices existing only when there is uniform background flow with the precise 
velocity to cancel the beta effect. 
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The formation of jets (localized, elongated energetic 
flows) in rotating 2D flow, such as planetary atmospheres 
and oceans, is a ubiquitous phenomenon. Although jet- 
like structures are often dictated by external forcing, or 
by boundary constraints, zonal (east-west) jets may also 
form via unforced evolution of essentially isotropic initial 
conditions in large systems where boundaries are unim- 
portant — the remarkable band structure of Jupiter is 
a famous example. Although the planetary rotation axis 
clearly distinguishes the zonal and meridional directions, 
its effect on the dynamics is rather subtle and there is 
no simple argument why this should lead to highly elon- 
gated structures. In particular, conserved quantities like 
energy E and enstrophy D,2, as well nonlinear advection 
effects, remain isotropic. 

Recently, a new, highly anistropic adiabatic invariant 
B, special to systems of interacting Rossby waves j^, has 
been argued to provide the required dynamical mecha- 
nism Q. In particular, for the ever larger scale flows 
generated by the well-known inverse cascade of energy, 
B increases strongly with scale if the energy spectrum 
remains isotropic. Conservation of B requires that the 
spectrum concentrate in the meridional direction, lead- 
ing to zonally concentrated flows. Well deflned Rossby 
waves exist only at mid-latitudes, so this offers only a 
partial explanation since jet formation is also seen in 
near-equatorial flows. 

In this work a very general, complementary approach 
is considered. Rather than trying to follow the very com- 
plicated turbulent dynamics from some given initial con- 
dition, one seeks to understand only the range of possible 
long-time steady states produced by free decay of such 
flows. In particular, equilibrium steady states are con- 
sidered, which are completely specifled by the values of 
all the conserved quantities. In standard 3D thermody- 
namic systems, often only energy and particle number are 
conserved. In the 2D flows of interest, there are an infi- 
nite number, providing a far greater range of interesting 
equilibria with nontrivial spatial structure 0, 0, S • 

In the following, an exact nonlinear PDE that deter- 
mines the equilibrium flows is derived, and solved in var- 
ious limits. It will be shown, quite generally, that the 
latitude dependence of the Coriolis force (the "beta ef- 



fect") destabilizes large coherent vortex structures: the 
only stable equilibrium flows are those that are zonally 
translation invariant, i.e., jet-like. Only in a background 
flow that precisely cancels the beta effect are coherent 
vortices stable. Interestingly, conservation of B, which 
requires a beta effect, also fails in this limit, and therefore 
does not hinder the formation of such structures. The 
equilibrium arguments are independent of all dynamical 
considerations, and hence do not depend on the existence 
of B, but the fact that they are so consistent suggests a 
relation. This will be a subject of future investigation. 

For generality, consider a class of equations (examples 
to follow) which may be written in the form 



9tQ -f V • VQ = 0. 



(1) 



The incompressible velocity field v = V x ^ = 
(dyijj, —dxijj). The stream function ^ is assumed related 
to the convectively conserved "charge field" Q through 
an energy functional Ti-IQ]: 



^(r) = SH/SQir). 



(2) 



An example is the Charny-Hasagawa-Mima (CHM) equa- 
tion in which Q = oj + kj^i/j + /, where lu = V x v = 
dxVy — dyVx — —V^ip is the vorticity. The Coriohs func- 
tion / = 2ilE sm[6{y)], where 9{y) is the latitude, and 
He = 27r/(24 hr) is the rotation frequency of the Earth, 
and coordinates are chosen so that the x-axis points east- 
ward and the y-axis northward. In this case. 



d\ / Sr'\Q{v) - /(r)]5(r,r')[Q(r') - /(r')]. 



(3) 

in which g is the Green function of the operator — V^-|-fc|j 
satisfying free slip and/or periodic boundary conditions. 
In the beta-plane approximation, f — fo + PfUi where 
/o = 217Esin(0o), /?/ = 2{nE I Re) cos{eQ), where Re 
is Earth's radius, and Oq is a reference latitude. The 
Rossby radius of deformation is i?o = l/^fi = c//, where 
c 3 m/s is the speed of internal gravity waves. In prin- 
ciple kfi = kii{y), but usually one sets — l/i?o(0). 
Only Pf then produces anisotropy. The free space 
Green frmction is the modified Bessel function, g(r, r') = 
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-fToC^filr — r'|)/27r, and linearizing Q produces the usual 
Rossby wave dispersion relation u = — (3 fk,j. / {k'j^ + fc^). 
The CHM equation (reducing to the Euler equation for 
kR = 0) is an approximation to the shallow water equa- 
tions in which propagating gravity waves are neglected. 
The surface height is then proportional to '0, and adia- 
batically follows the flow. This equation also describes 
drift wave plasmas, in the 2D plane perpendicular to an 
applied magnetic field, where if) is the electrostatic po- 
tential and uj the charge density. 

For later purposes, it is useful to define the Legendre 
transform of 7i, 



(4) 



in which is used to substitute ij} for Q. The relation 
may be inverted via the implied identity 



Q(r) - -5C/5^j{r). 
From one obtains the simple result 
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(5) 



(6) 



the first two terms of which may be recognized as the 
kinetic and potential energies. 

Since v • Vij: = 0, it follows that H is conserved 0. 
With periodic or free slip boundary conditions, it follows 
also from |^ that flh = J cPrh[Q{r)] is conserved for any 
ID function /i(cr), conveniently summarized by 



/ d\6[a-Q{v)], 



(7) 



conserved for any cr. One recovers Vlh = / d'7h{(j)g[a). 
Certain "momentum functionals" 



P = / d\X{r)Q{r) 



(8) 



are also conserved if Ti, has appropriate translation sym- 
metries . For rotational symmetry, the conserved (ver- 
tical component of) angular momentum corresponds to 
A — ir^. For translation symmetry along direction 1, the 

conserved linear momentum corresponds to A = 1 x r. 

For the CHM equation in the open beta-plane there is 
an additional adiabatically conserved quantity, which to 
quadratic order in takes the Fourier space form 
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d^k 
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6(k)Q(-k)Q(k), 



(9) 



in which Q(k) = + 

of - / = (-V2 + fc; 



fc|j,)7/;(k) is the Fourier transform 
The key properties are 0: 

(1) b{k) = 0[{kR/k)''] decays rapidly for k/ku > 1, and 

(2) for k/kfj <C 1, b(k) = 0(1) for k-space directions 
7r/3 < < 27r/3 (but vanishing for \0\ -> tt/2) while 
6(k) = 0{kji/k) otherwise. Property (1) implies that B, 



even more so than 7i, is a large scale quantity, insensitive 
to small scale variations in Q. Property (2) then implies 
that the inverse cascade must focus the small k part of 
the spectrum close to \9\ = tt/2, leading (through the 
curl relation with v) to large scale zonal flows ^ . 

When fcfl — > or /?/ ^ 0, the adiabatic conditions 
fail, and B is no longer conserved ^ , so this form should 
be used only at mid-latitudes. Since B is not exactly 
conserved, it will not be included explicitly in the equi- 
librium calculation, but its consistency with various pro- 
posed equilibrium states will be conflrmcd at the end. 

All statistical equilibrium information is contained in 
the free energy 



In / DQAg[Q]e 



-0{n[Q]^aoP[Q]) 



(10) 



in which (3 = 1/T is the inverse temperature, ap is a 
Lagrange multiplier for P, J DQ is a functional inte- 
gral over all possible conflgurations of Q, and Ag[(5] = 
ricr ~ 5cr[Q]) represents the infinite product of 
delta-functions required to impose the chosen values g{cr) 
of all the conserved integrals ^o-iQ] represented by the 
right hand side of lO. If Q{r) Qi is discretized on a 
grid with microscopic spacing a, then Liouville's theorem 
specifies the measure / DQ = lima^o Ilj I^oo 0- 

It transpires that may be computed exactly under 
the very general assumption that Ti. and P are insensi- 
tive to very small scale fluctuations in Q ^ H] - Due to 
fine-scale mixing, the equilibrium Q will fluctuate wildly 
from grid point to grid point, but its average Qo(r) over a 
small area P ^ (with I — > also at the end), will vary 
smoothly. If no(r, cr) is the equilibrium probability den- 
sity for finding a parcel of fluid with charge a in the area 
P about r. Then Qo{r) = J adano{r,a). It is assumed 
that H is smooth on scale I, and hence H[Q] — HlQo]- In 
this is provided by the smoothness of g: although it 
diverges at the origin, the logarithmic singularity is suf- 
ficiently weak that this assumption remains valid Q . By 
integrating out the small-scale fluctuations, which may 
be treated as independent from grid point to grid point, 
one obtains an entropy contribution. 



d^rdano{Y, a) ln[no(r, cr)]. 



(11) 



in terms of which the free energy is ^[tt-q] = TCyQa] — 
TS[nQ\. One now observes that if TS is to remain finite 
as a ^ (so that nontrivial equilibria are obtained in 
which entropy and energy compete) , one must adopt the 
scaling T = Ta^, /3 = with fixed /3 = l/T. 

To self consistently determine no(r, cr), Lagrange mul- 
tipliers are introduced via the Gibbs free energy 



d^r/x(cr)no(r, cr). 



(12) 



in which /i(cr) is used to tune g{<T). One may now freely 
extremalize Q over all no, constrained only by the nor- 
malization J dcrrio(r, cr) — 1, to obtain 



rio(i", cr) 



-^W(*o-a)„/3[/i(<T)-cr(*o-a)] 



(13) 
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where VtoCr) = ST-l/5Qo{r) is the equiHbrium stream 
function, a{r) = aQX{r), and 



W{t) = jln 



(14) 



Substituting H13|l back into l|12(l and using finaUy 
determines Q as a, functional of ^'q: 



(15) 



The latter is finally determined by the (minimum for T > 
0, maximum for T < 0) condition 



-SC/S-i'oir) = Qo(r) = F[*o(r) - a(r)], 



(16) 



where F{t) — —drW{T). Since W{t) is convex, i^(T) 
is monotonic. Using |SJ and the left hand side is 
Qq = (— V-^ + A:^)^'o + /, and (fT?)|l represents a nonhnear 
PDE for . The conserved quantities are set by the 
derivatives P = —dQ/dao — J cPrF{'i'o — a), g[a) — 
-6g/6fi{a)= Jd\noir,a). 

Substituting ifTCjl into (^3), one finds Q{r,t) = Qo{r + 
iaot) or Qo{r,0 + aot), depending on the choice of A. 
Fixed momentum solutions are not in general static, but 
require a background flow with velocity ao along the sym- 
metry direction. 

Equations H15|) , 1)16(1 are the fundamental results of this 
paper, providing a complete description of the fluid equi- 
libria for a rather general class of problems. Ultimately 
one is interested in coherent vortex formation, where Q 
is large in some compact region, and is much smaller out- 
side of it. One sees from (PJ that these are high energy 
configurations, and hence correspond to T < Q. 

For tractable applications "finite level systems" are 
useful: let g{a) = ^k=i ^kS{cr — ak), in which is 
the total area occupied by charge ak- One correspond- 
ingly requires only a finite number of Lagrange multi- 
pliers lik- e^^^-") = J2t=ie^'"'S{a - ak), and W{t) be- 
comes (the log of) a discrete sum. One obtains no(r, a) = 
ELi Pkir)S{a - CTfe), go(r) = J2k=i ^kPk{-r), where 



Pfe(i") = 7:^5 



(17) 



has spatial integral A^, and is therefore the equilibrium 
number density for charge Ok [To|- 

To illustrate properties of the solutions, consider the 
two level system, Uk — 0, ctq, in the beta-plane where a — 
aoy: p2 = (e^['^"(*o-°)-^l + 1)-^ is the Fermi function, 
where = /ii — /io is the chemical potential difference. 
Define the temperature scale To — ^'^o/^^r < 0' ^^'^ 
(^^0*0 - p)/2fo, qo = 2Qo/ao - 1, 
In terms of these scaled variables, 
the free energy and equilibrium equation take the form 



let t = T/To, po = 

P = {Px,Py) = kiiY 



r 
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FIG. 1: Scaled equilibrium charge density goq(f) (dashed line) 
and surface tension E(f) (solid line) versus reduced tempera- 
ture i = f /To. For f ^ 0, Qeq « 1 - 2e"^/*, E 1 - ■k'^P /12, 
while for t ^ 1, goq ~ \/3(l -t), E « [2(1 - t)f^'^. In- 
set: schematic coherent vortex with asymptotic charge den- 
sity a± = |(To(l ± gcq) inside and outside. Arrows indicate 
the direction of fluid flow, essentially parallel to the interface. 



^0 + ffoPy, with /lo = 1 - 2/o/(To + m/2To, 
The po-indepcndent eo(p) 



1^2 



90 



{-Vl + l)po 



tanh(po/i), 



(18) 



where h{p) 
90 

term is unimportant. The potential V{pq) — -^p^ 
t ln[2cosh(po/f)] is an even function with a single min- 
imum at po = for i > 1, and symmetric double minima 
at ±p°q(?) satisfying p^^ = i&-nk{pljt) = q°q for i<l. 

Equation (|18|) is a standard continuum model of a bi- 
nary fluid in a gravitational field go, composed of a mix- 
ture of heavy (go = 1, Qo = o")) and light {qq — Qo — 
0) particles. The |VpoP term is an attractive interaction 
between like particles which encourages phase separation 
at low temperatures, t < 1. For h 0, V{pq) — hpQ 
has a unique absolute minimum Pcq{t, h) — ~poq{t, —h). 
For t > 1, peQ is continuous through h — 0, while for 
i < 1 it jumps discontinuously between the heavy and 
light phases, ±p°q(f). This analogy makes it obvious that 
the stable equilibrium state must be vertically stratified, 
with po and qo increasing monotonically (lighter phases 
floating on denser phases) in the direction of go (south- 
ward if f3f = (3f + aok\ > 0). For i < 1 there will be 
a sharp interface centered at py — —ho /go (determined 
by p) of width Apy ~ t (see below) between the segre- 
gated phases. The analogy generalizes easily to multiple 
charge levels, which generate multicomponent fluid mix- 
tures with different combinations of separated and unsep- 
arated phases produced as i, pi are varied. However, go 
is unchanged, and in equilibrium the system must again 
be vertically stratified with phases density ordered along 
go- The argument also clearly generalizes to nonconstant 
I3f{y). A varying go{y) still produces vertical stratifica- 
tion. This establishes the claim that only equilibria with 
purely zonal flow are stable in the presence of a beta ef- 
fect, (3f ^ 0, and is completely consistent with inverse 
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scaled distance from interface 

FIG. 2: Scaled interface profiles qo{^,t), as a function of 
the scaled normal coordinate for a sequence of reduced 
temperatures from t — (discontinuous step, qo — sgn(^)) 
to t = 1 (completely flat, qo = 0) in steps of 0.1. For 
small i, qo ~ tanh(^/f) (hence width Ar « t/ka), while 
po ~ sgn(^)(l — e^'^'), remaining continuous even at t — 0. 
Near f = 1, go ~ Po ~ gcq tanh(^/25o) with diverging width 
Co = fcflAr = l/v'2(l -t). 

cascade arguments based on conservation of B. 

Only for /3/ = 0, hence ao = -0f/k% = 0(10 cm/s) 
does the gravity effect disappear In a background 

flow moving at speed ao, an effective isotropic /-plane is 
restored. Consider a large vortex region whose boundary 
is smooth on the scale of its width Ar{T) (see below), 
viewed here as a bubble of one phase inside the other 
(see inset to Fig. QJ. Coexistence requires T < Tq, and 
h = 0, i.e., /I = o-o(o'o — 2/o)/2A:|.. Near the interface, 
(|18|l becomes a ID equation in the coordinate £, = p ■ 
n normal to the interface, with gojPo i'Zeq [hence 



charges (y±{T) = \(Tq{1 ± g^q)], deep on either side. 

The free energy increment per unit length L of the 
interface, A^/L = {2\fo\^/'^lao)Y. yields the scaled the 
surface tension 

/oo 
d^id^Pof. (19) 
-oo 

In Fig. n numerical solutions for q^^, S are plotted (and 
exact asymptotics described), while in Fig. El scaled in- 
terface profiles qoi^) are plotted for several i. The true 
equilibrium solution is finally obtained by minimizing 
the vortex perimeter L at fixed area A (set by the to- 
tal charge Oi), yielding, not surprisingly, a circular vor- 
tex. These arguments again generalize to multiple charge 
levels, where one could in principle have more than two 
phases in simultaneous equilibrium, with vortices com- 
posed of a central core of one phase, ringed by one or 
more other phases. 

Since B is no longer conserved in the /-plane, there is 
no dynamical barrier to the formation of isotropic struc- 
tures, again establishing consistency between arguments 
based on the dynamics of the inverse cascade, and those 
based on purely equilibrium thermodynamics. It appears 
likely that conservation of i? is a microscopic reflection 
of the analogy to gravity-induced density stratiflcation. 
Since vortices drift across go, rather than accelerating 
along it, equilibration is necessarily far less direct than 
in physical binary fluid systems, and B may be one of 
the mechanisms controlling this. This idea will be inves- 
tigated more carefully in future work. 
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